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Abstract 



The presence of exciton phonon interactions is shown to play a key role in the exciton 
superfluidity. We apply the Landau criterion for an exciton-phonon condensate moving 
uniformly at zero temperature. It turns out that there are essentially two critical velocities 
in the theory. Within the range of these velocities the condensate can exist only as a bright 
soliton. The excitation spectrum and differential equations for the wave function of this 
^ ! condensate are derived. 

; PACS numbers: ll.35.Lk, 05.30.Jp, 63.20.Ls, 64.60.Ht 



The problem of critical velocities in the theory of superfluidity arose a long time ago when 
the experiments with the liquid He showed a substantial discrepancy with quantum-mechanical 



predictions. Later, the effect was analyzed and its phenomenological description was given (e.g. 



see H). The fact that the liquid He could not be treated as a weakly non-ideal Bose gas was 
believed to be the main reason for inconsistency of microscopic theory with experimental data. 

For a long time He has been the only substance where the superfluidity can be observed. The 
recent experiments with the dilute gas of excitons @, [FJ provide new possibilities for studying 



S-h ; different types of superfluidity. 

In this series of experiments the CU2O crystal was irradiated with laser light pulses of several ns 
duration. At low intensities of the laser beam (low concentration of excitons) the system revealed 
a typical diffusive behavior of exciton gas. Once the intensity of the beam exceeds some value, 
the majority of particles move together in the packet. Their common propagation velocity is close 
to the longitudinal sound velocity, and the packet evolves as a bright soliton. 

Some alternative explanations of the phenomena are known. One of them 0] implies that the 
bright soliton is a one-dimensional traveling wave which satisfies the Gross-Pitaevsky (nonlinear 
Schrodinger) equation 0] for the Bose-condensate wave function \P(x, £) 

ih?^- = -— + v^*m 2 (1) 
dt 2m K J 

with attractive potential of exciton-exciton interaction v < 0. 
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A quantitative treatment given in || provides an iterative solution for the Heisenberg equation 
with the use of perturbational methods. In this picture the second order interactions, neglected 
in the Bogoliubov approximation, contribute to the negative value of v. However, the influence of 
exciton-phonon interactions on the dynamics of the condensed excitons is not treated ||. 

Another interpretation is based on a classical model J!]] where the normal exciton gas is pushed 
towards the interior of a sample by the phonon wind emanating from the surface. Such an 
explanation has a discrepancy with the experiment, because the signal observed is one order of 
magnitude longer than the excitation pulse duration ||. 

In this Letter we give an alternative and, in our opinion, more intrinsic interpretation of these 
phenomena. We argue that it is a propagation of a superfluid exciton-phonon condensate which is 
observed experimentally. The presence of exciton-phonon interactions is crucial for a "soliton-like 
superfluidity". This interaction plays a key role when the propagation velocity approaches the 
longitudional sound velocity. 

We start with the Hamiltonian of the exciton-phonon system 

H = H ex + H ph + H int , H ex = ~ J #*(x)A^(x)dx+ ~ J ^*(x)^*(y) i/(x-y) ^(x)%)dxdy, 

Hph = I {^( x ) 2 + ^( V "(x)) 2 |dx, H iBt = J a(x-^*(x)^(x)(Vu(y))dxdy, (2) 

where ^ and u are the operators of the exciton and phonon fields correspondingly, c is the 
longitudional sound velocity and p denotes the mass density of the crystal. The field variables 
obey the following commutation relations 

#(x),**(y)] =M(x-y), fc(x), %(y)] = -i7%5(x - y), i,j = 1, 2, 3. 

In (U) we omit the terms with the transverse sound velocity, since the interaction of excitons with 
transverse sound waves is much weaker than with the longitudional ones. 

It is convenient to change the reference system when we consider a uniform motion of the Bose 
gas. The transition to the reference system moving uniformly with the velocity v = (v, 0, 0) is 
immediate. In new coordinates the classical field equations become: 

(ih^ + ^-A + ^f- Jv{* - y)|V(y, t)| 2 d2/ 3 ) t) = V(x, t) J a(x - y) (Vu(y, t)) dy (3) 



' d 2 d 2 2 d 2 

T^T - 2V— — + V 2 — j - C 2 



A) u(x,t) = ivy a(x- y )hKy,t)| 2 dy, (4) 



dt 2 dtdxi dx\ J ' p 

where ^(x, t) = ^{x\ + vt, x%, x^, i) exp (—imvxi/h). 

The l.h.s. of Equation @ is Galileian invariant, while the l.h.s. of (|4]) is Lorentz invariant. As 
a result, the system (H), (f|) is neither Galileian nor Lorentz invariant. As we will see later, it is 
due to this noninvariance that the effective potential of exciton-exciton interactions depends on 
velocity. 

Let us consider slowly varying solutions of the system (Q), (|j). In this (long wavelength) limit 
one can replace z/(x) and cr(x) by z/ <5(x) and Oo5(x), where u (> 0) and do denote the zero-mode 
Fourier components of the corresponding potentials. 
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Solving (^), one can express the bounded at infinity time- independent solution u(x) in terms 
of ^(x). The effective potential of the exciton-exciton interaction is obtained after substitut- 
ing this expression into (|3|). The phonon field makes this potential long-range, anisotropic and 
v -dependent. The potential becomes asymptotically attractive along the v-direction and asymp- 
totically repulsive in directions perpendicular to v. It follows that stability of the corresponding 
solutions ip = (f>(xi) exp(— iuot), U{ = 6nq(xi) is preserved under the one- dimensional reduction of 
the system (H), (|). The functions <j){x\), q{x\) obey the following equations 

/ h 2 d 2 \ ( a 2 \ mv 2 

X <f>[x x ) = i/ - , 9 ° 0(xi) 3 , A = -hu - -r- + Ca , (5) 



\2mdx 2 J \ (c 2 — v 2 )p / 

dq{xi) = c _ ^o0(^i) 2 /gx 
dx\ (c 2 — t> 2 )p' 

where the integration constant C is fixed by the condition q — > const as \x\\ — > oo. In the 
last equations assumed to be real. This choice does not change the result but simplifies our 
calculations. 

It follows from (|5|) that the effective potential becomes attractive when v exceeds the critical 
velocity 

v = J<?-(a%/v p). (7) 

But if v exceeds the sound velocity c, the potential becomes repulsive again. As for the solution 
varying in the direction n = {n\,n2,n^), ip — /( nx ) ex P( — iuot), u = ng(nx), the critical velocity 
is v (n) = vq/ cos(6 i ), | cos(0)| > v /c, where 9 is the angle between n and v. 

When v is less than the critical velocity (0), equations have the following stable sta- 

tionary solutions 

(i) = 0o = \Jn/V = const, u = const, and 

(ii) = o tanh (/30 o (xi - a)) , — jr— ^- = _ a °^ cosh" 2 (/3<f> (xi - a)) , (8) 

OXi (c z — v z )p 



(3 



mv \vl -v 2 \ ( \ 2 2 v 2 -v 2 a cj) 2 
A — T 7 ) T\ u o \ <Po — ^o9>o~5 o , ^ 



\ k \cr — v A \ \(c z —v z )p J c z — v z (c z — v z )p 



In (||,i) AT and V stand for the number of particles in the condensate and the volume of the system. 
When v exceeds Vq, we have only one stable stationary solution 

= 0o cosh" 1 (/?0 o (xi - a)) , d< i Xl = - , /"^ cosh" 2 (/30 Q (xi - a)) , 

axi (cr — v z )p 

-vS-'&tfA, c = a (9) 

To find the excitation spectrum of the system we expand the field operators near the proper 
classical solutions: 

$(x, t) = (0(xx) + x(x, t)) e- iu>ot , Ui(x, t) = 5 il q{x 1 ) + ^(x, t). 
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The Hamiltonian of the system can be written as follows 



H = H + hH 2 + (10) 

where H = H (<pe~ tuJot , q) stands for the classical part of H. It is important that H 2 is bilinear 
in x(x, t), fj(y:,t), whereas the linear terms are absent in ([H^) (since the classical fields satisfy the 
stationary equations (§),(§]))• From now on we are working in quasiclassical approximation and 
neglecting the terms of power greater than one (in h). 



The quasiclassical Hamiltonian ( |T0|) is reduced to the normal form 



H 2 = oJjblbj + const, 



hi, b* 
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bj 



0. 
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Indeed, since H 2 is a bilinear function of x, rj, the equations of motion are linear in field operators. 
They coincide with the corresponding classical equations (i.e. equations linearized around 

^(x,t) = 0(xi)exp(-iw o t), Ui(x,t) = Saqfa)): 



ih— H A 

at 2m 



X + Ca + 



<7n 



c 2 A - v x 



2 



v 2 )p 

d 2 



2u 



X - Vq<I>{x)\* 



+ 2v— — 



dt 2 



°"o, 



rj + ^Vi^ix + X*)) 
P 



a <j ) (x)(Vi 1 ) = 0, (12) 
0. (13) 



dx\ ' " dtdxi 

The quantities Ui in (|ll|) are characteristic frequencies of the system ([121) , flllf ). 

Let us consider the homogeneous Bose gas moving uniformly with velocity v < Vq. The 
condensate wave function is given by (§,i). The differential equations (0),(13) have constant 
coefficients so that the characteristic frequencies uj(k) are determined as roots of the following 
characteristic polynomial 



(ft 2 - c 2 k 2 ' 



h 2 (tt + vk 1 ) 



h 2 k 2 (K z k 



2m \ 2m 



+ 2*/ 



(c 2 — v 2 )p 



Ti 2 k 2 a 2 <p 2 k 2 
2m p 



0, (14) 



where ft = u;(k) — v k\ are the excitation frequencies in the crystal reference frame. In the limit 
cr — > one gets the Bogoliubov || spectrum huj[k) = ^ 4^ J^j^- + 2z/o0o) f° r the exciton gas 
as well as the free phonon spectrum ft = ck. When we switch on an exciton-phonon interac- 
tion, the spectrum u;(k) becomes v-dependent and the Landau criterion of superfluidity for the 
homogeneous Bose gas U has to be properly modified. 

The transition to the normal state occurs if there exists such k ^ for a given velocity that 
ft < 0, i.e. 

min(cj(k) — vk\) = 0. (15) 

k 

Analyzing (|T4D we obtain the value of the critical velocity vl for the homogeneous exciton-phonon 

gas 

v L 



y/Ml/m) (1 - (cx 2 /u pc 2 )) = Uu Q N/mV/c) v . 



The quantization near the translationally noninvariant classical solution (§,ii) in the region v < 
vq yields the same continuous spectrum u{k). The only new feature is that there appears a 
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bounded state at lu = in the v-direction. This fact has a simple explanation: the family of 
the solutions (|8|,ii) contains an arbitrary translation parameter a, which, in fact, is a collective 
coordinate. Differentiation of (|8|,ii) with respect to a gives then necessary time independent 
solution of fll^) )( |l3[) . This bounded state does not affect the quasiclassical excitation spectrum 
and contributes only to highest approximations (e.g. see ||). 

If the velocity v exceeds (0), the characteristic polynomial ( |14|) has complex roots and there is 
no stable constant solutions. The condensate (i.e. classical) wave function turns into the (bright) 
soliton (|9|) of the one-dimensional nonlinear Schrodinger equation (H). This solution decreases 
exponentially. This allows us to obtain the continuous spectrum from asymptotics of (|12D , (|13|) . 
We have 

hu(k) = A + - — 
2m 

for the exciton branch of the model, and 

oo(k) = ck + vki 

for the phonon branch. As in the previous case we get a bounded state at zero energy. We skip 
the question of existence of other bound states, since it is not essential for our purposes. 

The spectrum now has a gap in the exciton branch which is equal to A. In a sense, the situation 
is similar to the BCS theory: the exciton-phonon interaction makes the effective exciton-exciton 
potential attractive, and the excitation spectrum acquires a gap. 

The transition to the ballistic regime is accompanied by the symmetry breakdown: a new 
condensate wave function flj^) is no more translationally invariant. However, it contains a free 
translation parameter. We can interpret this as a phase transition of the second order. 

The value 0o is readily computed from the normalization condition J 0(x) 2 dx = N, and A is 
then obtained from @ 

mu^N 2 ( v 2 — v q 
8h 2 S 2 [c 2 -v 2 



In (fL6|) S denotes the packet cross-section in :r 2 :r3-plane. When v approaches the longitudional 
sound velocity c, the gap magnitude increases and soliton becomes more stable. The soliton energy 
can be estimated from ([|) 

E -^\24h 2 S 2 (c 2 -v 2 r [V +6VC +V °° bV ° V ) + 2 j + 

It follows from the last formula that E — > oo as v —>■ c. Roughly speaking, the soliton effective 
mass tends to infinity when its speed approaches the longitudional sound velocity. Then its motion 
is less subjected to the external forces. 

The onset of ballistical regime is determined by the condition v > vq, while the frictionless 
propagation of the packet is possible only if the condition min^ (o>(k) — vk{) > holds. The 
relation between the treshold speed v c for the frictionless propagation and the number N(v c ) of 
excitons in the packet at this speed follows from (|T5|) 

-,2 



c 2 - vj) _ u N(v c 

2 




v 2 -vl 2h S 
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where the function N(v)/S depends on the characteristics of the light source. As a consequence, 
v c > v and the velocity v c is also dependent on these characteristics. 

Thus the frictionless soliton motion is observed over the range v c < v < c. In principle, it is 
also possible to observe the onset of soliton propagation for some v close to v c with v < v < v c , 
but the packet should disintegrate at later stages of the propagation. 

It is easy to see that the solution @ is the most stable in the class of one-dimensional traveling 
waves moving uniformly with given v (> vq) and N. We argue that (P) is also the most stable 
solution in the class of all solutions with given v(> Vq) and N, because the effective exciton-exciton 
potential is attractive in v-direction and repulsive in the perpendicular directions. We would 
like to stress that effective one-dimensional solutions of three-dimensional nonlinear Schrodinger 
equations ([J) with attractive potentials do not have the similar properties. In particular, the 
stability of such solutions is doubtful |I0[ . 



When the crystal has no impurities, and the boundary friction is absent (the experiments with 
CU2O are carried out under these conditions, cf. 0,0), the propagation of solitons is observed 
at the critical velocity (|7|). 

The corresponding value of <7o is estimated from the deformation potential of CU2O. As a 
result, we obtain vq ~ 0.5 — 0.7c. These estimates are in agreement with the experimental data. 

To observe the critical velocity v c the presence of interior friction is important. Under these 
conditions, the relations between v c and the characteristic width of the soliton I is given by 

v c = Ti/ml. 

In the present work we have discussed the properties of the system at zero temperature. The 
extension of our results to finite temperatures seems to be a more difficult problem. 

We are very grateful to Y. Berest for valuable remarks and to L. Vinet and Y. Lepine for 
stimulating discussions. We are also grateful to A.Mysyrowicz for drawing our attention to the 
work j|] and useful comments. 
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